Inflation in loop quantum cosmology: 
Dynamics and spectrum of gravitational waves 
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Loop quantum cosmology provides an efficient framework to study the evolution of the Universe 
beyond the classical Big Bang paradigm. Because of holonomy corrections, the singularity is replaced 
by a "bounce." The dynamics of the background is investigated into the details, as a function of the 
parameters of the model. In particular, the conditions required for inflation to occur are carefully 
considered and are shown to be generically met. The propagation of gravitational waves is then 
investigated in this framework. By both numerical and analytical approaches, the primordial tensor 
power spectrum is computed for a wide range of parameters. Several interesting features could be 
observationally probed. 
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I. INTRODUCTION 



Loop quantum gravity (LQG) is a nonperturbative 
and background-independent quantization of general 
relativity. Based on a canonical approach, it uses 
Ashtckar variables, namely SU(2) valued connections 
and conjugate densitized triads. The quantization is 
obtained through holonomies of the connections and 
fluxes of the densitized triads (see, e.g., [l| for an intro- 
duction). Basically, loop quantum cosmology (LQC) is 
the symmetry reduced version of LQG (although it is 
fair to underline that the relations with the full theory 
are still to be investigated into the details). While 
predictions of LQC are very close to those of the old 
quantum geometrodynamics theory in the low curvature 
regime, there is a dramatic difference once the density 
approaches the Planck scale: the big bang is replaced by 
a big bounce due to huge repulsive quantum geometrical 
effects (see, e.g., Q for a review). Among the successes 
of LQC, one can cite: the excellent agreement between 
the trajectories obtained in the full quantum theory 
and the classical Friedman dynamics as far as the 
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density in much below the Planck scale, the resolution 
of past and future singularities, the "stability" of states 
which remain sharply peaked even after many cycles 
(in the k=l case) and the fact that initial conditions 
for inflation are somehow naturally met. The latter 
point is especially appealing as the inflationary scenario 
is currently the favored paradigm to describe the first 
stages of the evolution of the Universe (see, e.g., 0] for 
a recent review). Although still debated, it has received 
many experimental confirmations, including from the 
WM AP 7- Years results [3| , and solves most cosmological 
paradoxes. It is rather remarkable that, as will be 
explained in this paper, the canonical quantization of 
general relativity naturally leads to inflation without 
any fine tuning. Inflation would have been unavoidably 
predicted by LQC, independently of its usefulness in the 
cosmological paradigm. 

Two main quantum corrections are expected from the 
Hamiltonian of LQG when dealing with a semiclassical 
approach, as will be the case in this study mostly de- 
voted to potentially observable effects. The first one 
comes from the fact that loop quantization is based on 
holonomies, i.e. exponentials of the connection rather 
than direct connection components. The second one 
arises for inverse powers of the densitized triad, which 
when quantized become an operator with zero in its dis- 
crete spectrum thus lacking a direct inverse. As the sta- 
tus of "inverse volume" corrections is not clear due to 
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the fiducial volume cell dependence, this work focuses on 
the holonomy term only and derives, for the first time 
in a fully consistent way, the entire dynamics up to the 
explicit computation of the tensor power spectrum. The 
background evolution is first studied and a specific atten- 
tion is paid to the investigation of the inflationary stage 
following the bounce. Then, analytical formulas are given 
for the primordial tensor spectrum for either a pure de 
Sitter or a slow-roll inflation. Finally, numerical results 
are given for many values of the parameters of the model. 



II. BACKGROUND DYNAMICS 

In general, many different evolutionary scenarios are 
possible within the framework of LQC. However, all of 
them have a fundamental common feature, namely the 
cosmic bounce. As we will show, the implementation 
of a suitable matter content also generically leads to a 
phase of inflation. This phase is nearly mandatory in any 
meaningful cosmological scenario since our current un- 
derstanding of the growth of cosmic structures requires - 
among many other things- inflation in the early universe. 
It is therefore important to study the links between the 
inflationary paradigm and the LQC framework, as em- 
phasized, e.g., in [5|. 

The demonstration that a phase of superinflation can 
occur due to quantum gravity effects was one of the first 
great achievements of LQC [6| . This result was based on 
the so-called inverse volume corrections. It has however 
been understood that such corrections exhibit a fiducial 
cell dependence, making the physical meaning of the as- 
sociated results harder to understand. As reminded in 
the introduction, other corrections also arise in LQC, due 
to so-called holonomy terms, which do not depend on the 
fiducial cell volume. Those corrections lead to a dramatic 
modification of the Fricdmann equation which becomes 
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where p is the energy density, p c is the critical energy 
density, H is the Hubble parameter, and k = 8irG. In 
principle, p c can be viewed as a free parameter of theory. 
However, its value is usually determined thanks to the 
results of area quantization in LQG. Then, 
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where value 7 ~ 0.239 has been used, as obtained from 
the computation of the entropy of black holes Q ■ Should 
the inverse volume corrections be included, this would 
modify the background dynamics by some additional fac- 
tors. 

As it can easily be seen from Eq. ([1]), a general predic- 
tion associated with models including holonomy correc- 
tions is a bounce which occurs for p — p c . The appear- 
ance of this p 2 term with the correct negative sign is a 



highly nontrivial and appealing feature of this framework 
which shows that the repulsive quantum geometrical ef- 
fects become dominant in the Planck region. The very 
quantum nature of spacetime is capable of overwhelming 
the huge gravitational attraction. The dynamics of mod- 
els with holonomy corrections was studied in several arti- 
cles 0, lil-Hol] . In this paper we further perform a detailed 
and consistent study of a universe filled with a massive 
scalar field in this framework. The global dynamics of 
such models was firstly studied in Ref. Q. Recently, it 
was pointed out in Ref. @ that the "standard" slow-roll 
inflation is triggered by the preceding phase of quantum 
bounce. This general effect is due to the fact that the uni- 
verse undergoes contraction before the bounce, resulting 
in a negative value of the Hubble factor H. Since the 
equation governing the evolution of a scalar field in a 
Fricdmann-Robcrtson- Walker universe is 



2>H4> + m 2 4> = 0, 



(3) 



the negative value of H during the prebounce phase acts 
as an antifriction term leading to the amplification of 
the oscillations of field <p. In particular, when the scalar 
field is initially at the bottom of the potential well with 
some small nonvanishing derivative (j), then it is driven 
up the potential well as a result of the contraction of the 
universe. This situation is presented in Fig. Q] 
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1: Shark fin-type evolution of a scalar field for m = 
?7ipi. The (red) dot represents the point where the initial 



conditions in classical cosmology are usually set. 

To some extent, it is therefore reasonable to say that 
the LQC framework solves both the two main "problems" 
of the big bang theory: the singularity (which is regular- 
ized and replaced by a bounce) and the initial conditions 
for inflation (which are naturally set by the antifriction 
term). 

However, this shark fin evolution (see caption of Fig.[T]) 
is not the only possible one. In particular, a nearly 
symmetric evolution can also take place, as studied in 
Ref. [lj|. Those different scenarios can be distinguished 
by the fraction of kinetic energy at the bounce. When the 
energy density at the bounce is purely kinetic, the evo- 
lution of the field is symmetric. When a small fraction 



3 



of potential energy is introduced, which is the general 
case, the symmetry is broken and the field behaves as in 
the shark fin case. It is however important to underline 
that we consider only scenarios where the contribution 
from the potential is subdominant at the bounce, as it 
would otherwise be necessary to include quantum back- 
reaction effects [ll|. The effective dynamics would then 
be more complicated and could not be anymore described 
by Eq. ©. 

In order to perform qualitative studies of the dynamics 
of the model, it is useful to introduce the variables 



x :- 



m<fi 



and y := 



(4) 



Since the energy density of the field is constrained (p < 
p c ), the inequality 



y 2 < i 



(5) 



has to be fulfilled. The x 2 term corresponds to the po- 
tential part while the y 2 corresponds to the kinetic term. 
The case x 2 + y 2 = 1 corresponds to the bounce, when 
the energy density reaches its maximum. 

In Fig. [5J exemplary evolutionary paths in the x — y 
phase plane are shown. For all the presented cases, the 




FIG. 2: Exemplary phase trajectories of the scalar field with 
m — mpi. 

evolution begins at the origin (in the limit t — > -co), 
and then evolves (dashed line) to the point on the circle 
x 2 +y 2 = 1. Finally, the field moves back to the origin for 
t —> +oo (solid line). However, the shapes of the inter- 
mediate paths are different. The x = case corresponds 
to the symmetric evolution which was studied in Ref. 
(if the bounce is set at t = 0, the scale factor is an even 
function of time and the scalar field is an odd function). 
In this case, the field is at the bottom of the potential well 
exactly at the bounce (H = 0). This is however a very 



special choice of initial conditions. In the case x = V2/2, 
the potential term and kinetic term contribute equally at 
the bounce. In this case, both deflation and inflation oc- 
cur. However one observes differences in their duration. 
The third case, x = 1, corresponds to the domination of 
the potential part at the bounce. In this case, symmet- 
ric phases of deflation and inflation also occur (both the 
scale factor and the field being this time even functions). 
However in this situation, as well as in x = \/2/2 case, 
the effect of quantum backreaction should be taken into 
account. The dynamics can therefore significantly differ 
from the one computed with Eq. (fT|). 

In Fig. [3] we show some exemplary evolutions of the 
scalar field for different contributions from the potential 
part at the bounce. In Fig. 2J the corresponding evolu- 
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FIG. 3: Time evolution of the scalar field. Different evo- 
lutionary scenarios leading to a slow-roll inflation phase are 
displayed. The bottom (solid) line represents the symmetric 
case. The middle (dotted) line represents the shark fin-type 
evolution mostly investigated in this paper. The top (dashed) 
line corresponds to a larger fraction of potential energy. For 
all curves m = O.Olmpi. 

tions of the scale factor are displayed. It can easily be 
seen that the value of m ax increases with the fraction of 
potential energy at the bounce. Since the total energy 
density is constrained, </> ma x m ust satisfy 
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The values of </> ma x associated with different evolutionary 
scenarios were computed in 0, [i| [Io[ • The conclusion of 
those studies is that the necessary conditions for infla- 
tion are generically met. Only in the case of a symmet- 
ric evolution does the value of max become too small 
in some cases. In particular, for m = 10 _6 mpi one ob- 
tains </> m ax = 2.1mpi for a symmetric evolution. The 
corresponding number of e-folds can be computed with 

,2 

N ~ 271"-^-, which gives N ~ 28. By introducing a small 

m pi 

fraction of potential energy (as in the shark fin case), 
the number of c-folds can be appropriately increased. In 
addition to the usual arguments, this requirement is also 
set by the recent WMAP 7- Years results 0]. Based on 
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FIG. 4: Time evolution of the scale factor. Different evo- 
lutionary scenarios leading to a slow-roll inflation phase are 
displayed. The bottom (solid) line represents the symmet- 
ric case. The middle (dotted) line represents the shark fin 
type evolution. The top (dashed) line corresponds to a larger 
fraction of potential energy. For all curves m = O.Olmpi. 
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FIG. 5: The shark fin-type evolution of the scalar field for 
m = 10~ 3 7npi. Curves from bottom to top were computed 
for p c = 10~ 6 , 1CT 4 , 10~ 2 , 1, and 100 [m P1 ], respectively. 



III. GRAVITATIONAL WAVES IN LQC 



those observations, the value of the scalar spectral index 
was indeed measured to be ns = 0.963 ± 0.012. As, for a 
massive slow-roll inflation the relation 



1 m P1 
ns = 1 # 

7T 0^ 



(7) 



holds, one obtains </> bs = 2.9 ± 0.5mpi. Since the con- 
sistency relation <f> max > 4> t> s must be fulfilled, the sym- 
metric evolution with m = 10 _6 mpi (for which (/) max = 
2.1mpi < ^obs) is not favored by the WMAP 7- Years ob- 
servations. As already mentioned, higher values of cf> can 
be easily reached if some contribution from the poten- 
tial term is introduced (this supports the shark fin sce- 
nario). The number of e- folds will therefore be naturally 
increased in this way. However it remains bounded by 

above: since N ~ 27r-^j-, Eq. ((6]) leads to the constrain: 



N < 



(8) 
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The value of the parameter p c can be fixed by Eq. j2]). 
However, this expression is based on the computation of 
the area gap as performed in LQG. This, in general, can 
be questioned [l2|. In particular, in the framework of 
reduced phase space quantization of LQC, the value of 
p c remains a free parameter [l3j]. Moreover, a particu- 
lar value of the Barbero-Immirzi parameter (imposed by 
black hole entropy considerations) has been used. There- 
fore, the value of p c can, in general, differ and it is worth 
investigating how the variation of p c can alter the dynam- 
ics of the model. In particular, we have studied how the 
shark fin scenario can be modified by different choices of 
p c . In Fig. [5l the evolution of the field is displayed as a 
function of the value of the critical energy density. As 
expected, the larger p c , the higher the maximum value 
reached by the field. It can be seen that <^ max approaches 
the usually required value ~ 3mpi for p c ~ m P1 , making 
the whole scenario quite natural. 



Although quite a lot of work has already been devoted 
to gravitational waves in LQC [3, this study aims at 
treating, for the first time, the problem in a fully self- 
consistent way with an explicit emphasis on the investi- 
gation of the spectrum that can be used as an input to 
study possible experimental signatures. 

The equation for tensor modes in LQC is given (see, 
e.g., M) by 



— K + 2aH—K - V 2 K + mlK = 0, 



drj 



dr) 



(9) 



where h l a are gravitational perturbations, 7/ is the confor- 
mal time and the factor due to the holonomy corrections 
is given by 



UnGa 
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This factor acts as an effective mass term. For conve- 
nience we introduce the variable 



aha 



ah 6 



(11) 



where h\ = — h% — /i®, h\ 
the Fourier transform 

u(x,r,) = 



h\~h^. Then, performing 



d 3 k 



(2tt) 3 

one can rewrite the equation as 



Uk(n)e 



where k 



2 



d 2 

^2 U k(T?) + [fc 2 + m 2 ff ] U k (7?) = 0, 



k • k and 



'eff 



= a 2 



P-- 3 P 
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(13) 



(14) 
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It is worth underlining that the final expression of m e g 
has no explicit dependence upon the critical energy den- 
sity p c . In Eq. 
However since 



l4|) . both TOq and a /a depend on p c 
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(15) 



the factors depending on p c cancel out precisely. This is 
perhaps not a coincidence and this could exhibit the con- 
servation of classical symmetries while introducing the 
quantum corrections. 

The next step consists in quantizing the Fourier modes 
Uk{v)- This follows the standard canonical procedure. 
Promoting this quantity to be an operator, one performs 
the decomposition 

Mv) = fk(v)k + mv)bl k , (16) 

where fk(v) is the so-called mode function which satis- 
fies the same equation as Uk(??), namely Eq. ([T3| . The 
creation (&£) and annihilation (Sk) operators fulfill the 
commutation relation [b^, b^] = <5( 3 )(k — q). 

The problem is now shifted to the resolution of a 
Schrodingcr-like Eq. ([T^| which can be used to compute 
the observationally relevant quantities. In particular, the 
correlation function for tensor modes is given by 



\h%(x,v)h b a (y,v)\0) = 



dk„ ,, . sin kr , . 

T ?MM)— d7) 



where Vt is the tensor power spectrum and |0) is the 
vacuum state. In our case, Vt can be written as 



7MM) 



64ttG k 3 
a 2 (r]) 2^ 



\Mv)\' 



(18) 



This spectrum is the fundamental observable associated 
with gravitational wave production. As will be shown in 
the next sections, very substantial deviations from the 
usual shape are to be expected within the LQC frame- 
work. 



Basically, the phcnomcnological parameters entering the 
model are therefore: 

• Vi — the beginning of the inflation. 

• At? — the width of the bounce. 

• fco — which is approximately equal to the value of 
m e ff at the bounce (when H = 0). It can therefore 
be related with the energy scale of the bounce. 



• v — which is related to slow-roll parameter e by 
v = J I + 3e = | + e + 0(e 2 ), where e <C 1. 



For the considered model, we have /cq > 0. This comes 
from the fact that we consider the particular shark fin- 
type of evolution where the bounce is dominated by the 
kinetic energy term. Therefore when i/>i [see Eq. (U)] , 
Eq. (fT4| simplifies to m 2 cS = a 2 K0 2 /6 > 0, leading to 

fcp Si Tilsit = ^bounce) ^ 0. 

A matching should be performed between the three 
considered phases. It can be done, as displayed in Fig. HI 
with transition matrices defined as follows: 



M 



fk(v) /:(»?) 



(20) 



where the Wronskian condition implies 

W(/ fe ( 7? ),/ fe *(r ? )):=detM = i. (21) 
The inverse of the transition matrix is then given by: 



M 



(22) 



ra 2 eff 



IV. ANALYTICAL INVESTIGATION OF THE 
POWER SPECTRUM 

In this section we perform analytical studies of gravita- 
tional wave creation in the scenario previously described. 
In particular, we derive approximate formulas for the ten- 
sor power spectrum at the end of inflation. In the next 
section we will compare this result with numerical com- 
putations. 

In the considered model, the evolution is split into 
three parts: contraction, bounce and slow-roll inflation. 
For this model, the effective mass square is defined as 
follows 



err 



(V) 







for 
for 
for 



77 < 77j - A??. 

v > m- 
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k 2 



2_ 

2 

'II 







n 



FIG. 6: Evolution of the effective mass used in the analytical 
approximation [Eq. (|19[) ], On this plot, e is set to zero as 
an example. The dashed line represents the case without a 
bounce. The points where the transfer matrices are computed 
in our model are also indicated. 
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The three first transition matrices are: 



M, 



M, 



where 



V2k 

-MJfe-A.?) 



V2k 



/2n 



/2f2 

n p <fl(>j,-A ? ) 



2<> 



/2I2 



n = \/fc 2 + fcg. 



> (23) 
,(24) 
(25) 

(26) 



In the last region, mode functions can be written as 
fk(v) = a k9k (v) (27) 

where 

ffkfa) = V^JleW'+WHWi-kri), (28) 

H v (x) being a Hankel function of the first kind. The 
mode functions g k (rj) correspond to another decomposi- 
tion of the field u^rj) in the form: 



(29) 



The creation (d k ) and annihilation (dk) operators fulfill 
the commutation relation [dk, ojj = 5^(k — q). Because 
decompositions ([T6")) and (|2T)|) arc equivalent, based on 
Eq. (|27p and on the Wronskian conditions for the mode 
functions f k and one obtains: 



(30) 



which corresponds to a Bogoliubov transformation with 
coefficients a k and f3 k . Because of the commutation re- 
lation of the creation and annihilation operators we have 
|c*ft| 2 - \Pk\ 2 = 1. It is clear from Eq. §U§ that if 
Pk =/= particles are created from the vacuum, just be- 
cause 6k 1 0) = /^dlklO)- By matching the three regions, 
the unknown coefficients a k and [3 k can be determined: 



" 6k " 




a k pi 




dk 













a k 



= M 4 _1 M 3 M 2 ^ 1 M 1 



(31) 



Ma 



e ifc(Ari ~"- ) (Ocos[Ai)Q]-tfcsin[A»)ri]) 
\/2kQ 

e ,fc(A„-r,,)(_ tfc cos [A77n]-f2sin[A77n]) 
V2k 



where M4 is given by 
M 4 = 



9k(v) 9k(v) 
d v 9k(v) dv9* k (v) 



(32) 



the mode function gk being given by Eq. (|28|) . In the 
special case corresponding to a dc Sitter inflation (e = 
and v — |), the mode functions given by Eq. (|28|) 
simplify to the Bunch-Davies vacuum 



— ikr\ 

V2k 



1 



1 

kr] 



(33) 



In general, the amplitude of the mode function during 
inflation can be written as 



\fk 



\g k \ 2 \a k 



fa? 



m(al/3 k g* k )$tg k 



(34) 



As we are interested in the spectrum at the end of infla- 
tion (77 — >■ 0~), the approximation 



H 



(35) 



holds and, based on this, one can easily see that for a 
slow-roll inflation (e <C 1): 



lim 



O(e). 



(36) 



Therefore, the leading order contribution from Eq. (|34[) 
becomes 



lim \f k \ 2 



\g k \ 2 \a k 



P k 



(37) 



With this approximation, the tensor power spectrum at 
the end of inflation takes the form 



r T (k) 



16 ( H_\ 2 ( k 



7r \mp\ / V aH 



-2c 



\a k -p k f 



(38) 



The coefficients a k and k arc computed from Eq. (|31[) . 
Since the resulting expression for \a k -- (3 k \ 2 is very long, it 
is not explicitly given here. It exhibits the correct ultra- 
violet (UV) behavior, namely limfe_ i . 00 \a k — f3 k \ 2 = 1. 
Therefore, the UV spectrum simplifies to 



V T (k -> 00) 



16 ( H V f k 



(39) 



V=Vi 



In Fig. [71 spectra, as obtained from Eq. (|3"8"]l . are dis- 
played for different values of ko and normalized to the 
usual non-LQC corrected spectrum. In Fig. [8l the width 
of the bounce A77 is varied. In both cases, e is vanishing. 

The main features that can be drawn from those plots 
are the following: 

• The power is suppressed in the infra-red (IR) 
regime. This is a characteristic feature associated 
with the bounce. 

• The UV behavior agrees with the standard general 
relativistic picture. 

• Damped oscillations are superimposed with the 
spectrum around the "transition" momentum fc„ 
between the suppressed regime and the standard 
regime. 
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and of the background: 




FIG. 7: Analytical tensor power spectra, normalized to the 
non-LQC-corrected spectrum, for three different values of ko 
in the e = case. The parameters are: ko = (solid line), 
ko — 1.5 (dashed line), ko = 3 (dotted line), r\i = — f, and 
A77 = 1. 
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FIG. 8: Analytical tensor power spectra, normalized to the 
non-LQC-corrected spectrum, for two different values of Arj 
in the e — case. The parameters are: Arj = (solid line), 
Arj = 10 (dashed line), ko = f , and r\i = —I. 



• The first oscillation behaves like a "bump" that can 
substantially exceed the UV asymptotic value. 



• The parameter ko basically controls the amplitude 
of the oscillations whereas Arj controls their fre- 
quency. 



NUMERICAL INVESTIGATION OF THE 
POWER SPECTRUM 



To perform a more detailed analysis, we have also 
fully numerically solved the system of coupled differential 
equations which leads to both the evolution of the modes 



where 



d 2 h 
dt 2 
dH 
~dt 
da 

~dl 

d<f> 

~dt 

dn^ 

~~dT 



P = 



dt 

\ k (p + p) 
Ha, 



a z 6 



2^-1 



p) 



2a 6 



m 
~2 



and p = 



2a 6 



m 
T 



(40) 

(41) 
(42) 
(43) 
(44) 

(45) 



are respectively the energy density and pressure of the 
scalar field whereas 7r^ is the momentum. 

To compute the evolution of the modes, the initial con- 
dition was assumed to be the Minkowski vacuum 



fk — 



-ikrf 



2k 



(46) 



This approximation is valid for the subhorizontal modes. 
Therefore, in the numerical computations we have 
evolved only modes that were subhorizontal at the ini- 
tial time. 

In Fig. [9j the analytical spectrum Eq. (|38|) evaluated 
as explained in the previous section is compared with 
the full numerical computation. The overall agreement 
is very good with slight deviations due to subtle dynam- 
ical effects. The UV tilt associated with the slow-roll 
parameter is perfectly recovered. The values of parame- 
ters H, ko and e were determined from the evolution of 
the background. In turn, the parameters rji and A77 were 
fixed to fit the numerical data. 

The mass of the scalar field is, of course, the key phys- 
ical parameter of this model. The canonically chosen 
value (around 10 _6 mpi) may not be especially meaning- 
ful in this approach as the standard requirements of in- 
flation are substantially modified by the specific history 
of the Hubble radius. This value is nonetheless still the 
mostly preferred one. 

In Fig. [TQ1 the spectra computed for three different 
mass values are displayed. As expected, the UV value of 
the spectrum scales as m 2 , since during inflation Vt ~ 



H 2 



It is also clear that the region of oscillations 



becomes broader while lowering the value of m. 

In Fig. 11 11 we show how the spectrum is modified by 
different choices of p c . It is clear that increasing p c leads 
to an amplification of the spectrum. The dependence is 
however not very strong. As it was shown in Section ILT1 
the increase of p c leads to an increase of the field displace- 
ment (/> m ax- This dependence was shown to be rather 



weak. Since T T ~ H 2 



,2i2 



, the increase of cj> due to 



the dependence upon p c will result in an amplification 
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0.001 



FIG. 9: Comparison of numerical and analytical spectra 
Eq. (|38[) for m = 10~ 2 mpi. In the IR region the spectra 
behave as Vt oc k 2 while in the UV region they behave as 
Vt oc k~ 2e , where e <C 1 is the slow-roll parameter. Here: 
H = 0.037mpi,e = 0.0246, k = 0.037m P i,7/; = -750, and 
Ar/ = 10. 



m=10 m P i 




that increasing p c amplifies the oscillatory structure. 

The numerical investigations performed for this work 
have shown that the quantity R defined as 



R 



basically evolves as 



V T (k = K) 

'inslanda 
/ T 



{k — kit ) 



V m J 



0.64 



(47) 



(48) 



where £;* is the position of the highest peak in the 
power spectrum and 7-r tandard (fc) is a standard inflation- 
ary power spectrum [see e.g. Eq. ([39]) ] which overlaps 
with VT(k) for k — y oc . The function (H51 was ob- 
tained by fitting the numerical data in the mass range 
m = 5 • 10~ 5 TOpi ... 10 _1 mpi. Because of numerical in- 
stabilities, it was not possible to perform computations 
for lower values of the inflaton mass. The numerically 
obtained values of R together with the approximation 
given by Eq. (j48]l are given in Fig. [T2l This parametriza- 



R 



1000 



10- 



icr 4 o.ooi o.oi 



FIG. 10: Numerically computed power spectra for m = 
10 -4 , 10~ 3 , 10~ 2 ?npi (from bottom to top in the UV range) 
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FIG. 11: Numerically computed power spectra for p c = 
10~ 4 ,10~ 2 ,10 mpi (from bottom to top in the UV range) 
with m — 10~ 3 mpi. 



of the power spectrum. This is in agreement with the 
numerical results. From Fig. II H it can also be noticed 



0.0001 



0.001 



0.01 



0.1 



m 



FIG. 12: Ratio defined by Eq. (|47[) as a function of inflaton 
mass in Planck units. Dots are values obtained from the nu- 
merical computations. The straight line is the fit given by 
Eq. {48]). 

tion is useful for phenomcnological purposes. Interest- 
ingly, R can become very high for low values of the mass 
of the field. This partially compensates for the lower 
overall normalization of the spectrum and can become 
a very specific feature of the model. In particular, for 
the mass m ~ 10 _6 mpi (which is the value preferred by 
some estimations), extrapolating the relation (|48)) leads 
to R w 8000. If the relation still holds in this range, 
the effect is very significant, and could have important 
observational consequences. 

Finally, to make basic studies easier, we performed a 
rough parametrization of the full spectrum: 



16 / H 
Vt = — [ 

7T V ?77pi 



■ k i 



-2< 



1 + iK/kf 



AR-2 



1 + {k/Kf 



(49) 
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leading to 



V$ s 



16 / H 

ir \mpi 



1 



1 + (K/kf 



1 



4i?-2 



1 



(k/K) 

(50) 

in the specific case of de Sitter inflation. In both cases, 
the classical behavior is recovered in the limit k — > oo. 
The point for introducing the R factor the way it was 
done becomes clear when calculating the value of the 
spectra at k = fc*. For a modified de Sitter spectrum 
[Eq. fl50j)], we get 



Vf(k = K) =R^(JL 

TT \TOpi 



(51) 



Thanks to the relation (|48|) . the number of the free pa- 
rameters can be decreased in a phenomenological analy- 
sis. 

As shown on Fig. 1131 this formula correctly reproduces 
the main features, namely the IR power suppression, the 
bump and the UV limit. Oscillations are missed but due 
to momentum integration there is little hope that they 
can observationally be seen on a cosmological microwave 
background (CMB) spectrum. 




1CT 5 1CT 4 0.001 0.01 
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FIG. 13: Comparison of the numerical spectrum for m — 



10 mpi with formulas (|49|) and (|50|) . The solid (blue) line 
corresponds to 
to (|5U)). 



49l) while the dashed (red) line corresponds 



To conclude this section, we have schematically repre- 
sented the evolution of the Hubble radius (i?H := 
together with the physical modes, in Fig. Q31 This helps 
to understand the shape of the obtained spectra. 

We consider the modes that are initially (at time t\) 
shorter than the Hubble radius. For those modes, the 
normalized solution is given by the Minkowski vacuum 
fk = e-' lk ^/V2k. Therefore, the initial power spectrum 
takes the form Vt ~ k 3 \fk\ 2 ~ k 2 . Starting from the 
largest scales, the modes cross the Hubble radius. This 
is possible since the Hubble radius undergoes contraction 
faster than any particular length scale. While crossing 
the horizon, the shape of the spectrum becomes frozen 
in the initial Vt ~ k 2 form. Then, the modes evolve 



Hubble 
radius 




length 




ti contraction f 2 £3 inflation £ 4 



time 
— f>- 



FIG. 14: Schematic picture of evolution of the Hubble radius 
(solid line) and of the different length scales (dashed lines) 
for the considered model of the universe. Different times are 
distinguished: t\— time when the initial conditions are set; 
ti— bounce (H = 0); ts — beginning of inflation; ti— end of 
inflation; t$— present epoch of dark energy domination. 



through the bounce (at time t%) until the beginning of 
inflation (at time £3). The main consequence of the tran- 
sition of modes through the bounce is the appearance of 
additional oscillations in the spectrum. This issue was 
studied in detail in Ref. 16[ , where the spectrum at time 
£3 was calculated for the symmetric bounce model. Af- 
ter the bounce, modes with wavelengths shorter than A* 
start to reenter the Hubble radius. The superhorizon 
modes A > A* (k < fc*) hold the k 2 spectrum, with how- 
ever some oscillatory features due to the bounce. Modes 
with A < A* (k > k*) cross the horizon again during 
the phase of inflation. For them, the spectrum agrees 
with the standard slow-roll inflation spectrum Vt ~ k~ 2t 
where e <C 1. The small tilt is due to a slow increase of 
the Hubble radius. Contributions from different modes 
are then slightly different. At the end of inflation (at time 
£4) the spectrum is therefore suppressed (Vt ~ k 2 ) for 
k < fc* and exhibits the inflationary shape (Vt ~ k~ 2e ) 
for k > k*. The spectrum is also modified by the os- 
cillations due to the bounce. This corresponds to the 
computations of this paper. The particular mode with 
wavelength A* (wave number k*) should be studied in 
more detail. The size of this mode overlaps with the 
size of the Hubble radius at the beginning of inflation: 
fe* ~ a(t3)H(t^)- The physical length A* at the scale fac- 
tor a(t) is therefore equal to A*(i) ~ a(t)/[a(t 3 )H(t 3 )]. 
This scale grows with the cosmic expansion and it is 
crucial, from the observational point of view, to deter- 
mine its present size (at time t§). The case drawn in 
Fig. [2] corresponds to a present size of A* greater than 
the size of the horizon. This is indeed rather unlikely 
that the present value of A* is below the size of horizon 
just because the spectrum of scalar perturbations should 
then exhibit deviations from the nearly scale invariant 
inflationary prediction. Up to now, there is no observa- 
tional evidence for such deviations. A remaining possi- 
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bility would however be that the (slight) observed lack of 
power in the CMB spectrum of anisotropics could be due 
to the effects of the bounce. However, the present size 
of A» would then be comparable with the size of horizon. 
This leads to the question: why should those two scales 
overlap right now? This is rather unnatural, and would 
lead to a new coincidence problem. However, as it was 
estimated in Ref. |{|, these two scales can indeed over- 
lap in the standard inflationary scenario for quite natural 
values of the parameters. There is therefore a glimpse of 
hope that the scale A* is at least not to much bigger than 
the size of horizon. This could allow us to see some UV 
features due to the bounce as the oscillations also affect 
sightly the inflationary part of the spectrum. These are 
however secondary effects and it is not clear whether they 
were not smoothed away during the radiation domination 
era. Moreover, in the region where those effects could be 
expected, errors due to the cosmic variance become sig- 
nificant. This is an unavoidable observational limitation 
which cannot be bettered, even by the improvement of 
resolution of the future CMB experiments. 

Another limitation in studying the effects of LQC 
comes from the fact that the derived modifications can 
also appear in other bouncing cosmologies. In particu- 
lar, within the model of quintom bounce, the discussed 
effects of suppression and oscillations were also pointed 
out [Il,[l9j]. The amplitude of tensor perturbations at the 
peak was however not predicted to be as high as in LQC. 
An additional amplification on the very large scales was 
also predicted in the quintom model. Despite these dif- 
ferences, at the observationally accessible low scales, the 
effects due to the LQC bounce and the quintom bounce 
are mostly indistinguishable. Therefore, complementary 
observational methods have to be proposed to distinguish 
between such models. A possible distinction could be 
given e.g. from the analysis of non-Gaussianity produc- 
tion within LQC. 

VI. CONCLUSIONS 

This study establishes the full background dynamics in 
bouncing models with holonomy corrections. Although 
this has already been claimed before, we confirm that 
due to the sudden change of sign of the Hubble param- 
eter, inflation is nearly unavoidable. In this paper, we 
have considered a particular model of inflation where the 
content of the universe is dominated by a massive scalar 
field. We have investigated the details by both analyti- 
cal and numerical studies the primordial power spectrum 
of gravitational waves. It exhibits several characteristic 
features, namely a Vt °c k 2 IR power suppression, os- 
cillations, and a bump at k*. In the UV regime, the 
standard inflationary spectrum Vt oc k~ 2e is recovered. 
The primordial tensor power spectrum transforms into 
B-type CMB polarization. The performed investigations 



therefore open the window for observational tests of the 
model, in particular through the amplification which oc- 
curs while approaching k — > fc* . The observed structures 
correspond to the UV region in the spectrum. If the 
present scale A* ~ 1/fc* is not much larger than the size 
of horizon, then the effects of the bounce should be, in 
principle, observable. In particular, one should expect 
amplification, rather than suppression of the B-type po- 
larization spectrum at the low multipolcs. The suppres- 
sion for k < fc* becomes dominant at the much larger 
scale, probably far above the horizon. While the B-type 
polarization has not been detected yet, there are huge 
efforts in this direction. Experiments such as PLANCK 
[T3|, BICEP [H or QUIET [U are (partly) devoted to 
the search of the B mode. Even with present observa- 
tional constraints, one can already exclude some evolu- 
tionary scenarios and possible values of the parameters, 
in particular the inflaton mass m and position of the 
bump fc* in the spectrum. We address this interesting 
issue elsewhere [22( ■ There are also still several points to 
study around this model: 

• How is the scenario modified when quantum back- 
reaction is taken into account (in particular when 
the potential energy of the field in not negligible at 
the bounce)? 

• How is the power spectrum modified by inverse- 
volume terms in this framework? Although the 
background dynamics should not be fundamentally 
altered, the spectrum could be significantly modi- 
fied. 

• How do those results compare with models deal- 
ing with classical bounces (see, e.g., (23[)? If the 
IR power suppression is probably a generic fea- 
ture of bounces, the detailed features are model- 
dependent. 

Together with the known success of LQC (The singu- 
larity resolution, the correct low-energy behavior, etc.), 
the facts that 1) inflation naturally occurs and 2) obser- 
vational features can be expected from the model, are 
strong cases for loop cosmology. Those two points are 
the main results of this paper. 
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